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Abstract 



'^ I Spreading of different types of fluid on substrates under an impressed 

Q ■ force is an interesting problem. Here we study spreading of four fluids, 

O . having different hydrophilicity and viscosity on two substrates - glass 

and perspex, under an external force. The area of contact of fluid 
and solid is video-photographed and its increase with time is mea- 
^^ . sured. The results for different external forces can be scaled onto a 

f^ I common curve. We try to explain the nature of this curve on the basis 

CN ' of existing theoretical treatment where either the no-slip condition is 

used or slip between fluid and substrate is introduced. We find that 
of the eight cases under study, in five cases quantitative agreement is 



in 

Q I obtained using a slip coefficient. 
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1 Introduction 

Wetting and spreading [H [21 [3] are very common everyday phenomena, but 
there are still a lot of open problems in this field. A simple example is - 
what happens when a fluid is forced to spread under an externally impressed 



force ? There is no fully satisfactory answer to be found in the otherwise 
comprehensive literature on the topic, though such experiments are routinely 
used to characterize viscoelastic materials. 

In the present paper, we report a set of very simple experiments on the 
so-called 'squeeze film' and try to find a theoretical framework to explain 
our results. The problem is briefly described as follows - a small droplet of 
a liquid, placed on a solid plate will take the shape of a section of a sphere. 
Under equilibrium, the well known Young-Laplace relation holds, involving 
an equilibrium angle of contact 6'e.The dynamic condition while the fluid is 
spreading, i.e. the area of contact of the fluid and solid is increasing in time, 
is also well studied. There are a number of excellent reviews summarising 
studies on - the rate of spreading, the forces responsible and other interesting 
features [H El IH [3] . We now look at a simple extension of the problem - what 
happens if the fluid is forced to spread beyond its equilibrium area, under a 
constant weight placed on top of it. A prehminary report was published by 
our group in [6]. In this paper we study four different fluids on two different 
solids, so as to have different combinations of viscosities of fluids as well as 
different degrees of hydrophilicity of the fluids and solids. The results are 
analyzed on the basis of viscous energy dissipation with and without slip [Till]. 

We find that while the no-slip result gives a fair quantitative estimate of 
the contact area as function of time for only one case, it is possible to get 
a better agreement in 5 of the 8 cases, on introducing a slip coefficient. Of 
course, this involves introduction of a free parameter (the slip coefficient) 
whereas the no-slip analysis did not. The shp coefficient appears to depend 
more on the solid surface than the fluid, and is an order of magnitude lower 
for the perspex plate. The results for spreading of the same fluid on different 
surfaces is quite different. There are 3 cases where the agreement between 
theory and experiment is unsatisfactory and cannot be improved by intro- 
ducing slip, so other effects are also at play. We describe the experimental 
setup followed by the theory and present detailed results in the third sec- 
tion. Finally we compare experimental results with theory and discuss the 
implications. 

2 Experimental 

The experimental fluids are castor oil, linseed oil, glycerol and ethylene glycol. 
The choice of fluids is motivated by the following idea. While castor oil and 
linseed oil are insoluble in water, with a low dielectric constant e the other two 
fluids ethylene glycol and glycerol having high dielectric constants, are soluble 
in water, (e = 4.67, 3.35, 47.70 and 43.0 respectively for castor oil,linseed 



oil, glycol and glycerol). Again, one fluid belonging to each type according 
to hydrophilicity, has high viscosity(/i), namely castor oil and glycerol(/i for 
castor oil=986-451 cp for temp 20-30c,glycerol is 1490-629 cp for temp 20- 
30c), while the others have much lower viscosity(/i for glycol is 19.9-9.13 cp 
for temp 20-40c , linseed oil=33. 1-17.6 cp for temp 20-50c). 

A constant volume Vf of the fluid is measured with a micropipette and 
placed on the lower plate. The upper plate is then placed on the drop and a 
load of M kg. (M = 1 to 5) is placed on it. The pair of plates are ~ 1 cm thick, 
one pair is made of float glass and the other of perspex. A video camera, 
placed below the lower glass plate records the slow increase in area. The video 
clips are analyzed using Image Pro-plus software. We find that the area of the 
drop in contact with the plate increases with time at first. Then after some 
time it reaches a saturation value and does not increase significantly with 
time any more. The time of saturation and area of saturation are different 
for different fluid-solid combinations. We plotted the area vs. time graph 
for all fluid-sohd combinations. Fig(IT]) and fig([2]) show respectively area vs. 
time graphs for Castor oil and glycol respectively on perspex for loads from 
1-5 kg. Interestingly when the data for area are divided by the load mass 
M, the curves come very close to each other, almost collapsing to a single 
master curve. Again, for castor oil on perspex area increases for a long time. 
It takes several minutes to saturate approximately and continue to increase 
very slowly even after 5 minutes. For glycol on both the substrates, on the 
other hand, the area increases very fast and saturates in a few seconds. Video 
photos are taken initially, at the rate of 10 frames/sec, in this case. 

3 Theory 

We discuss in this section standard formulations for this kind of problem. The 
'squeeze flow' usually discussed in chemical engineering literature is different 
from our experiment, ours is a constant volume squeeze, where whole mass of 
fluid taken is always confined within the two plates. In the usual experiment 
a constant area setup is used, where both plates are immersed in the fluid, 
the distance between the plates as function of time is measured as a force is 
applied on the upper plate. 

First we outline the steps worked out for the no-slip case, i.e. assuming 
that the layer of fluid next to the static plate has zero velocity. The analysis 
roughly follows [7]. 



3.1 No-slip 

We assume the fluid to be Newtonian and incompressible with constant vol- 
ume Vf. Ours is a constant force set up, the force Mg acts on the fluid layer. 
The velocity first builds up, reaches a maximum and then starts to reduce as 
it squeezes the liquid film. The details are presented in [B], here we consider 
the velocity decreasing regime only. We equate the loss of potential energy 
of the loaded plate with the work done against the viscous resistance. With 
h and r representing the instantaneous thickness and radius of the film re- 
spectively, p the pressure and pa the atmospheric pressure. In moving over 
a distance dh the work done is given by 



dh (p — pA)2Trrdr (1) 

The loss of potential energy is 

Mg{h -{h- dh)) = Mgdh (2) 

Now the pressure p at any point r from the center of the upper plate for 
any given distance h, where the velocity of the upper plate V = —dh/dt can 
be calculated as 

Therefore, we can write, 

/ Mgdh = dh {p- pA)2iTrdr (4) 

Substituting for (p — pa), we get, 

Mg= ^ ^ ^-2TTrdr (5) 

Which , on integration gives, 

3/if ? dh , , 

^ 27r/i5 dt ^ ' 

3u(t;f)^ 11, 



The right hand expression depends on the impressed force Mg through the 
area of contact A{t). If we divide both sides by M, the LHS is independent 
of M. The RHS is now a function of A{M), M and t. We call this function 
G{M, t). We calculate the RHS from experimental data and compare with the 
LHS. If the principles on which the derivation is based are correct, the RHS 
calculated from experimental results should be independent of M. That is 
the curves G{M, t) vs. t should collapse to a single curve and be numerically 
equal gt. 

so, we test the result 

gt = G{M,t) (10) 

by plotting the LHS and RHS separately on the same graph and check 
how well they agree. Of course the agreement can be expected only for the 
time interval before the area of contact saturates, since the LHS goes on 
increasing linearly. 

3.2 Analysis of Experimental Results assuming no-slip 

The results for glycerol on glass are shown in figQ. We see that for this 
combination assuming no slip gives an acceptable quantitative agreement 
between the Ihs and rhs of eq fllOp . This calculation involves no adjustable 
parameter. It may further be noted that the data for different loads collapse 
almost to a single line, on scaling by M. However, for the other fluid-solid 
combinations, such quantitative agreement is not found. We try to see now 
if this discrepancy may be attributed to the effect of shp at the fluid and 
solid interface. 

We follow the equation given by fS] listed in P] for this situation. 

3.3 Effect of slip 

Relative motion between a solid surface and material contacting it (wall 
slip) is a phenomenon observed with many materials [1]. In the last section 
we assumed that the layer of the liquid touching the lower fixed plate has 
zero velocity. So that the applied force is entirely carried away by viscous 
dissipation. But experiments show that there may be cases where the layer of 
the liquid having contact with the substrate may have finite velocities, which 
may affect the spreading of the liquid drop. In these cases, shear stresses at 



the sample-plate interfaces may vanish (for perfect slip) or reduce (for partial 
slip). The amount of slip depends not only on the liquid tested, but also on 
the material and roughness of the substrates. We introduce partial slip in 
the analysis as follows. 

3.4 Partial slip 

No slip and perfect slip are two limiting cases of slip phenomenon. To explain 
our results, we take partial slip into consideration. Stefan [8j gave a slip law 
of the form , 

V=I3*T (11) 

where, /3 is the slip coefficient and r is the wall shear stress . The value 
of /3 is zero for no slip and is infinity for perfect slip cases. It may have 
any value in between, for partial slip cases, depending on the liquid and the 
substrate. 

With the introduction of /3,the force equation becomes, 

3/it;? f dh 
Mat = — -^ / — 14 

Dividing by M. 

gt = R{M,t) (15) 

where, the integral can be evaluated as 

R{M, t) = -1/k^ * [kyh^ + 3ky2h^ - ln{l + k/h)]f^ (16) 

with k = 6/5 /x 

The above force equation, on integration, gives us the function R{M,t) 
similar to G{M, t) for the no-slip case. We fit the slip coefficient /3 so as to 
get the best fit between the left and right sides of equation [151 the results 
are shown in the next section. 



3.5 Results with partial slip 

Table(Il]) shows values of /3 for different fluid-solid combinations. We find 
that the values of (3 ranges from 10"'^ to 10~^ for glass and from 10~^ to 10~® 
for perspex. So glass has a larger /3 than perspex. For a definite substrate 
the values of (3 are different for different liquids. In Table 1 we also show 
the initial contact area Aq (before loading) and the viscosities and dielectric 
constants of the fluids and substrates. 

Figures (jl]) and ([5]) show respectively the graphs comparing R{M, t) and 
gt for ethylene glycol on glass and for castor oil on perspex. The coefficients 
P for the best fit are given in table 1. 

While five of the eight curves can be fit quite well with a finite /3, in the 
remaining three cases a discrepancy remains. Here the fit becomes worse for 
any finite (3. This is shown for glycol on perspex in figure®, the other two 
cases are glycol on glass and linseed oil on perspex. 

4 Discussion 

Out of eight cases studied, we find that the theoretical framework considered 
here can account for five cases fairly well with one free parameter, the slip 
coefficient (3. This is not bad considering the complexity of the problem and 
the simplifications introduced. 

Scrutiny of the table 1, reveals several interesting features. Let us look 
at Aq, the area of contact with only the plate on the fluid drop and no 
extra weights, we consider this as the reference condition. The mass of the 
perspex plate is 250 gm and the glass plate is about 700 gms. It is seen 
from spreading vs. load data (see [6]) that extrapolating to zero load does 
not change this reference point significantly. We see that for glass Aq ranges 
from 2-5 cm^, while for perspex we have values less than 1 and the highest 
is 1.4 cm^. So equilibrium spreading on glass is higher in all cases. However, 
under load fluids spread to a much greater extent on perspex, compared to 
glass. Looking at the percentage (of Aq) spreading for the 5 kg. load, we 
see that for glass it is always less than 50 %, whatever the fluid, while for 
perspex the area increases by 100 to 400 %. So, where the initial spreading 
is large, the effect of load is less and vice versa. 

The highest spreading with load is observed for castor oil on perspex, 
moreover castor oil shows a creep like behavior and continues to spread for 
very long times, though the rate slows down. Possibly visco-elastic effects 
also play important role and should be taken into consideration. 

The three cases which cannot be explained by the present theory are all 



for the low viscosity fluids ethylene glycol and linseed oil. This treatment 
focussing on the viscous effect works better in case of high viscosity. We have 
tabulated the dielectric constants along with other data in table 1, because 
earlier work [9] showed that the polarizabilties of the fluid and substrate are 
important in spreading phenomena. But here we find no significant correla- 
tion except that the maximum spreading (under load) is for the low dielectric 
fluid on the low dielectric substrate. 

To conclude, spreading of different types of fluids on different solids is a 
very complex and interesting problem, still open for exploration. Very simple 
experimental apparatus provides a lot of food for thought in this area. The 
spreading and flow characteristics are routinely used for classifying soft mate- 
rials in food technology and chemical engineering, so a deeper understanding 
will be useful. 

Further experiments and theoretical study taking into account effects of 
rheology, surface roughness and ambient conditions is required 
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Table 1: The table shows properties of the fiuids and substrates used, /i 
represents the viscosity, e^ and ei, the solid and liquid dielectric constants, 
S the surface tension. Results from the experiment are - the initial area of 
contact Ao, the percentile spreading of area under loads of 1 and 5 kg., the 
slip coefficient /3 which gives a best fit for the spreading master curve. 
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Figure 1: This figure shows area vs. time graph for Castor oil on perspex. 
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Figure 2: This figure shows area vs. time graph for glycol on perspex. 
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Figure 3: This figure tests gt=G(M,t) results for glycerol on glass 
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Figure 4: This figure tests gt = R(M,t) results for glycol on glass. Introduc- 
tion of slip coefficient /3 gives a good agreement. 
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Figure 5: This figure tests gt = R(M,t) results for castor oil on perspex. 
Introduction of slip coefficient /3 gives a good agreement. 
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Figure 6: This figure shows that gt does not equal R(M,t) for glycol on 
perspex, for any positive value of f3. 
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